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Визначено поняття та особливості АР гравітаційного моделювання. Обробку його даних 
на коротких профілях запропоновано проводити за співвідношенням зі швидкоспадними 
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методикою Двуліта. 
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К обоснованию мониторинга. На фоне 
очевидного глобального технологического про- 
гресса в условиях зкономического спада на 
количество, методику и возможности осуще- 
ствления крупномасштабньх геофизических 
работ влияют две тенденции. С одной сторо- 
нь, возросшая точность и производительность 
гравиметрических сьемок (зтот метод при при- 
емлемой точности остаєтся недорогим поиско- 
во-разведочньм решениєем в общем комплек- 
се геофизических исследований вследствие со- 
вершенствования аппаратурьт и топогеодези- 
ческого обеспечения -- систем СРЗ5). С другой 
стороньк, заметно резкое снижение обьема из- 
мерений. 

Первое повлекло за собой необходимость 
пересмотра методов обработки полученньх 
данньхх, в частности более точньг6й учет по- 
правок Буге |Бьгчков, 2007). Последнее на фо- 
не необходимости обнаружения более глубин- 
ньх источников аномалий! влечет за собой 


| Возможности методов регуляризации в реше- 
ний задач построения сложньх разрезов при со- 
временном уровне модельньсх представлений о гео- 
логической среде близки к технологическому пре- 
делу. 
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пересмотр методики измерений с целью уче- 
та тонких особенностей гравитационньх ано- 
малий "малой кровью", т.е. без усложнения ма- 
тематического аппарата, методик измерений и 
увеличения материально-технических затрат. 
Зти особенности можно " подсечь", если за- 
действовать дополнительную переменную -- 
время. 

В зтом плане мировьеіе "трендьт" геофизи- 
ческих наблюдений постепенно тяготеют к не- 
прерьвному 4Р мониторингу |Вабіачіа еї аї., 
2008) исследуемой площади, изучая зволюцию 
гравитационного поля за время зксплуатации 
участка либо на протяжений интервала его рез- 
кой динамической активизации. Следуєт под- 
черкнуть, что в западной литературе | Сеорбу- 
зіс5, 2008 | под термином "вариации сильг тя- 
жести" имеют в виду временную разность ре- 
альньх аномалий на ограниченньжх площадях, 
источниками которьх вьступают обьекть! с 
бьстроизменяєемой глубинной динамикой, в то 
время как в отечественной литературе за зтим 
понятием зарезервированьт слабьге квазипери- 
одические флуктуации на сверхдлинньх про- 
филях, пересекающих области контрастньх со- 
временньх вертикальньгх движений блоков зем- 
ной корь. 
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Сама идея мониторинга уже витала в воз- 
духе |Лоссовский, 2006)? и некоторье ее при- 
ложения в определенной степени иллюстри- 
рует работа | Юргин, 2006). Тем не менее, в 
указанньх источниках не упоминаєтся о не- 
обходимой, по нашему мнению, подоплеке мо- 
ниторинга -- вариациях сильг тяжести в том 
пониманий, которое в них заложено работа- 
ми Г.Т. Собакаря, и их зависимости от ряда 
малойнтенсивньгх природньгх и техногенньх 
факторов. 

Применение повторньїх наблюдений в не- 
которой области в гравиметриий для нас срав- 
нительно неново (его злементь  применялись 
в созданиий региональньгх плотностньх моде- 
лей Кавказа, о чем упоминаєтся в работе |Але- 
ксидзе, 1985|, атакже |Юркина, 1978|), хотя 
непрерьвньми во времени их назвать трудно, 
в то время как в сейсмометрий зто давно об- 
щепринято. 

Отметим, что при организации любь(х вьі- 
сокоточньх измерений, в которьігх сила тяже- 
сти -- отдельньгй параметр, ее изменение сле- 
дует обязательно исследовать и учитьгвать |Со- 
бакарь, 1972). 

О вариациях сильг тяжести. В работе 
ІСобакарь, 1972| впервьге в отечественной ли- 
тературе проанализированьт различнье со- 
ставляющие вариативной части гравитацион- 
ного поля и вьщеленьт неприливньг6е квазипе- 
риодические вариации (КПВ ) сильг тяжести 
(самплитудой в пределах трех-пятикратной по- 
грешности измерений). Для их измерений бь- 
ла создана новая методологическая | неодно- 
родная модель Земли как в физико-химичес- 
ком, так и в знергетическом плане) , методи- 
ческая ( КПВ максимальньг при пересечений 
разновозрастньх тектонических структур и 
в областях контрастньх современньтх верти- 
кальньхх движений) и метрологическая (наблю- 
дения на сети оптимальной плотности и кон- 
фигурации несколькими приборами с вьісоким 
козффициентом надежности и сравнениєе гра- 
диентов КПВВ и наблюденньх аномалий) база. 

Сами КПВ, очевидно, тесно связаньт с зн- 
догенньми процессами формирования и раз- 
вития плотностньжх неоднородностей земной 
корьт и мантиий, хотя не исчерпьваются ими. 


ЗС.82-- "если современньже физико-матема- 
тические теории строят и исследуют свом функ- 
циональнье связи и последующие смьісловьіе рядь 
на территории преимущественно пространства, 
то зксперимент зондируєт пока еще не подвласт- 
ное теории время". 
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В ходе перераспределения вещества внутри Зем- 
ли меняєтся ее момент инерций, а как следст- 
вие -- ротационньй режим и сила тяжести. 
При зтом перераспределениий нарушаєтся фи- 
гура равновесия Земли, что изменяєт силу тя- 
жести. В ходе изостатического восстановления 
фигурьг равновесия сила тяжести снова меня- 
ется. Зта гамма процессов и генерирует КПВ 
силь тяжести Земли. Обратимая часть процес- 
сов создает периодическую часть вариаций, 
анеобратимая часть -  непериодическую часть 
вариаций, которая и формируєт стационарное 
гравитационное поле Земли (гравианомалим). 

М.В. Лубковьгм |Лубков, 2011| отмечена кор- 
реляция крупномасштабньсх плотностьньтх не- 
однородностей мантии и приливньх парамет- 
ров Земли, характеризирующих ее приливное 
деформированиєе. Сами приливнье параметрьі 
(так назьтвваємьге числа Лява и Шида) находят- 
ся на основе расчета релаксационньгх ампли- 
туд гравипотенциала на поверхности Земли 
и ее ядра. В связи с зтим возможно наличие 
кросскорреляции КГІВ и упомянутьгх динами- 
ческих параметров Земли, что требуєт отдель- 
ного изучения. 

Невзирая на слабую интенсивность, анома- 
лий КПВ сильттяжести вьшявляются вполне уве- 
ренно вследствие особенностей поведения вре- 
менной кривой КГІВ и тесной обратно пропор- 
циональной корреляции с кривьіми вертикаль- 
ньх движений земной корьх. Последнюю счи- 
тают | Собакарь, 1972| следствием общности 
процессов в верхней мантий, влияющих на 
КПВ и вертикальнье движения. Суммарная ве- 
личина КГІВ является суммой суперпозиций ва- 
риаций различного происхождения, знака, пе- 
риода и амплитудьт. Зта суммарная амплиту- 
да в конечном итоге и определяет зволюцию 
гравитационного поля Земли, вьізванную зво- 
люцией неоднородностей ее корьг и верхней 
мантий. 

Ввиду сказанного, вариации сильт тяжести 
следуєт учитьгвать в интерпретации результа- 
тов АР гравиметрии -- сцелью введения соост- 
ветствующих поправок в гравитационнье сьем- 
ки разньх лет, в долговременную прецизион- 
ную топосьемку, уточнения реологиий исследу- 
емого участка и др. 

Основа мониторинга. Определим основ- 
нье понятия. Назовем гравитационньм мони- 
торингом ряд периодически повторяємкеіх в 
реальном времени непрерьівньх на протяже- 
ний фиксированного временного отрезка (рис. 
1) микрогравиметрических измерений грави- 
тационного поля и обработку его данньх с уче- 
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Рис. 1. 11-дневньй ряд измерений силь тяжести: 
а-- без фильтрации, б -- с полосовой фильтра- 
цией и в -- с поправкой за приливньй зффект. 


том влияния окружающей средьтг и сферьт при- 
менения. 

Величина временного отрезка зависит от ка- 
чества измерений, мерьі неопределенности ре- 
зультатов наблюдений, динамики (амплитудьт 
и частоть) поля. 

Физической основой гравиметрического мо- 
ниторинга являєтся непрерьгвная связь дина- 
мики гравитационного поля и параметров сре- 
дьг - - ундуляциям сантиметрового диапазона 
соответствуют вариации силь тяжести в не- 
сколько мГал. Если деформация дневного рель- 
ефа определенного участка есть прямим след- 
ствием приповерхностного распределения масс, 
гравитационньй мониторинг можно применять 
для исследования разуплотнения и флюидно- 
го режима зтого участка. 

Пространственное распределение вариаций 
вертикальной производной У, потенциала си- 
ль тяжести прямо коррелирует с площадньтм 
распределением плотностей, а временньте ва- 
риации У, четко определяют вертикальнь!е ва- 
риации насьшщения флюийдов. 

Аппаратная основа мониторинга средьт пред- 
ставляет собой совместньве крупномасштабньте 
измерения превьшений рельефа по данньтм СР5 
и абсолютньгх значений силь тяжести (сотни 
стационарньгх пунктов на сотни км2). Значи- 
тельно дешевле обходится гравиметрия с от- 
носительньгкм учетом силь тяжести, но ей при- 
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сущи серьезнье ограничения -- привязка к 
опорной сети и необходимость синхронного 
учета "сползания нуля". Тем не менее, в на- 
ших условиях именно такие измерения при 
надлежащем методическом обеспечений име- 
ют наибольшую перспективу. На некоторьжх 
участках возможно применять непрерьтвньте 
измерения силь тяжести в скважинах, преоб- 
ладающие в разрешающей способности назем- 
ную сьемку благодаря большей близости к ис- 
точникам и устранению приповерхностньїкх вли- 
яний. Для компенсации влияния температурь 
и других внешних факторов используют неко- 
торье наработки из аппаратурной базьт мор- 
ской гравиметриий. 

О вьщделений слабого сигнала. В целом, 
в исследованиях, посвященньх АР гравимони- 
торингу | Сеорпубісз5, 2008), отмечена амплиту- 
да сигнала в диапазоне от 20 мГал |М/іШата5- 
ЗДопе5 еї а!., 2008| до 80 мГал |Сеіпдб еї аї., 
2008). В границах действующих вулканов ам- 
плитуда сигнала возрастаєт до 300--600 мГал, 
ав пределах водньгх бассейнов - 200--250 мГал 
ГОаміз еї а!., 2008) при нелинейном учете "спол- 
зания нуля". Вследствие повторяемости изме- 
рений, фиксирующей невязку значений сильт 
тяжести в точке наблюдений, можно оценить 
точность определения полезного сигнала на вре- 
менном отрезке. 

С целью вьделения слабого сигнала на фо- 
не помех в гравиметриий применяют методьт 
корреляционного анализа и вьїтисление не- 
которьхх составляющих градиентов по прямьтм 
измерениям силь тяжести?. Авторьт |Вайадіа 
еї а!., 2008| вькделяют полезньг6й сигнал путем 
расчета разностной аномалий Дд д - Ай п 7 
-УП, - Ада 7 Ду, где А - - разность меж- 
ду двумя соседними отсчетами во времени, 
у є- - 308,6 мГал / м -- поправка в свободном 
воздухе, П, - - вертикальное смещение (отно- 
сительное поднятие- -опусканиє), Л9 цер - - ано- 
малия Буге деформации (вклад от изменения 
обьема из-за сжатия средні вокруг источника, 
что подразумевает смещение плотностньх гра- 
ниц в неоднородной срелде); Лд у з 2ЛСруфб, з 
- 42965, -- влияние грунтовьхх вод (действиєе 
гидротермального фактора). Кроме того, ис- 
пользуется привязка к независимьтм наблюде- 


ЗТочно так же, как в морской гравиметричес- 
кой сьемке применяют вьічисление вертикальньтх 
и горизонтальньх градиентов силь тяжести и их 
комплексную интерпретацию. Зта методика интер- 
претации опробована на морских нефтегазовьхтх 
структурах ПЮргин, 2006). 
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ниям на контрольном пункте, расположенном 
вблизи исследуемого участка. 

Определенньге улучшения слабого сигнала 
на фоне помех можно получить и путем ста- 
тистической фильтрации влияния температур- 
ньїх флуктуаций с помощью алгоритмов с не- 
четкой логикой |Апао, Сагропе, 2004). Если не 
хватаєт длиньг ряда наблюдений, используют 
полиномиальную аппроксимацию низкого по- 
рядка |Вопуаїої еї а!., 1998| с соответствующей 
"калибровкой" порядка полинома. Кроме того, 
используют сравнение результата фильтрации 
с данньми, полученньтми с близкорасположен- 
ного контрольного пункта от участка наблю - 
дений. 

Примечательно, что как и в случае измере- 
ний КПВ сильг6 тяжести, многие авторьі |Сео- 
рбузвісз, 2008| отмечают корреляцию слабого 
сигнала с вертикальньми подвижками и, пред- 
полагая разночастотность напряженности поля 
сильт тяжести и вариаций помех, применяют 
для усиления полезного сигнала частотную фильт- 
рацию. Аналогию можно усмотреть опять же в 
морской гравиметриий. На равнинньх участках 
сигнал имеет мальгїй градиент, и такой подход, 
по нашему мнению, малозффективен. В таком 
случае необходимо учитьвать влияние помех 
иньжми методами -- изменением геометрии се- 
ти наблюдений, временной фильтрацией, вьі- 
числением производньх ит. п. |Собакарь, 1972). 

О методике мониторинга. Вьїіскажем не- 
которье соображения об организации мони- 
торинга. Его методическую сторону, видимо, 
целесообразно основьівать на достижениях |Со- 
бакарь, 1972), а изучение метрологических ню- 
ансов -- на примере геодинамического поли- 
гона с готовой инфраструктурой и базой гео- 
физических измерений (например, одного из 
использованньх в исследованиях |Двуліт, 1999|). 
Так, для мониторинга условного месторожде- 
ния необходима сеть детальной рядовой сьем- 
ки масштаба 1 : 10 000 наровном полигоне 10 х 
х 10 км? по стандартной методике с миниму- 
мом опорньгх пунктов. Согласно работе |Гра- 
виразведка ..., 1990, с. 106), 21 профиль (по 20 
точек с шагом 50 м) с контролем 30 9; от об- 
щего количества рядовьгх точек (свьше 13 тьс. 
пунктов измерений). Организация полигона и 
наблюдения на указанной сети за месяц обой- 
дутся приблизительно в 100 тьс. грн. без уче- 
та стоимости и амортизации аппаратурьі и дру- 
гих накладньх расходов. 

Недостаток инфраструктурьг и оборудования 
увеличит затратьт на организацию мониторин- 
гового полигона на неопределенную стоимость. 
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Использование цифровой регистрации позво- 
лит за счет упрощения схемьі сьемок обойти не- 
которье трудности в начальньх зтапах сьем- 
ки и ускорить создание цифровой модели обь- 
екта. Вьшход из финансовотехнологического ту- 
пика -- сотрудничество разньх институтов с 
совместньтм использованием инструментария, 
персонала, методик и общим планом исследо- 
ваний на общем обт-екте (зтим путем интегри- 
рованного мониторинга (соПарогайує дгаміїу пеї- 
угогК) идут многие западнье компаний). 

Следуєет отметить, что применениєе класси- 
ческой схемьт измерений на регулярной сети 
пунктов и последующий перерасчет значений 
по известному интегралу Пуассона пригодно 
для региональньх исследований |Раміз еї аї., 
2008|, но в локальньгх условиях, где чаще все- 
го применяют гравимониторинг | Болотнова, 
20071, имеет ряд недостатков |Дубовенко, 2002. 
К тому же, иногда по ряду причин, организа- 
ция регулярной сети невозможна, а пересчет 
из нерегулярной сети на регулярную являєт- 
ся задачей более сложного порядка, чем об- 
ратное восстановление строения средбі по дан- 
ньм измерений. 

Решение обратньхх задач гравиметриий с дан- 
ньшми в псевдорегулярньгх сетях с примене- 
нием моделей средьт типа "бесконечньгй про- 
филь" приводит к появлению плохо обуслов - 
ленньгх систем линейньгх уравнений, генери- 
рующих бессодержательньге результатьїі. Вви- 
ду зтого и исходя из преимущественно крат- 
кой длиньт реальньж профилей, целесообра- 
зен альтернативньігй подход. 

Для интерпретации измерений на коротких 
профилях предложена система линейньх ин- 
тегральньїхх уравнений с бьістро убьюввающими 
ядрами |Дубовенко, 20021: 


"м 00 2(х) о з-за 


х | 84 (2) | сов ог денна (9), 
ДБА 20и (х) 

2 1 

5пз1(Х) 2 У(х) - 65 


со -1 


5-5 (Хх), 
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Со (0) 2 50 (9 як 5с (х), побу о, 


С учетом зтого, методика, изложенная в ра- 
боте |Болотнова, 2007), остаєтся действенной 
в определенньгх условиях (региональньгй фон 
-- полином 1-й степени; известньтт плотнос- 
ти и положенця границ тяготеющих тел на 
поверхности; зти тела близки или имеют об- 
щие контакть). Она, в частности, предпола- 
гает построение пространственной плотност- 
ной модели средьт в три зтапа: разделение ано- 
малий от тяготеющих тел и определение зф- 
фективньгх глубин их залегания и квазиплот- 
ности (нулевое приближениє), определение ис- 
тинньх глубин и плотностей тел по решению 
плоской обратной задачи (первое приближениєе) 
и окончательное решение в среде АРС-3Р. 

Автором предлагаєтся в методике |Болот- 
нова, 2007| использовать программньй комп- 
лекс |Старостенко и др., 2004), разработанньтй 
на основе работьт | Старостенко, Легостаєва, 1998|, 
и подборку программ, полученную в исследо- 
ваний |Дубовенко, 2002. 

К интерпретации данньх. Целью мони- 
торинга является оценка глубиньі к источни- 
ку аномалий и изменения обьема по данньтм 
деформации рельефа, которьге требуют знания 
приповерхностного распределения масс (по 
данньгкм гравиметриий). Деформации земной по- 
верхности получают по данньтм спутниковой 
геодезий СР5, имеющих ряд преимуществ пе- 
ред традиционньми геодезическими метода - 
ми?. Приповерхностнье неоднородности стро- 
ения средьт (карстьк, плавунь(к, зоньг обводне- 
ния и разуплотнения), сложная структура пло- 
щади (складчатость, солянокупольная тектони- 
ка, разломіі), факторь поглощения полезного 
сигнала (температура, инструментальнь!е вли- 
яния) ограничивают зффективность мониторин- 
га, не снижая его практического значения. 

В работе |Ваначііа еї а!., 2008; Раміз еї аї., 
2008) не учтена особенность вариаций грави- 
тационного поля: колебания величиньї его про- 
изводньх зависят от флуктуаций малоийнтен- 
сивньїх геофизических факторов (аномальньх 
атмосферньх масс (рис. 2), снеговьіх масс, 
уровня грунтовьгх вод, лесистости и измене- 
ний рельефа (рис. 3) вследствие техногенной 


й Независимость от времени суток, погодьі, авто- 
матика, непрерьівность, полнота, надежная привяз- 
каксети. 
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Рис. 2. Изменение гравиполя Западной Украийньт 
вследствие перемещения атмосферньхх масс |Дву- 
літ, 1999. 


деятельности |Двуліт, 1999). Учитьгвать зти зф- 
фекть можно путем внесения соответствую- 
щих поправок при решений прямих задач гра- 
виметрии на исследуемьх территориях (пред- 
полагаєтся, что вследствие длительного мони- 
торинга участка приблизительно известно его 
строениєе). 

К примеру, влияние притяжения атмосфер- 
ньжх масс в центральной сферической зоне 
(сферическая трапеция 1? х 1?) вокруг пункта 
наблюдений на величину гравитационного по- 
ля определяет формула | Двуліт, 19991 540 - 
- 0,269АР , где АР -- аномалия атмосферно- 





23) 
о 
о « 
о 


Рис. 3. Изменения гравиполя Западной Украйньт 
вследствие перемещения масс в недрах |Двуліт, 
19991. 
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го давления; влияние снеговьгїх масс в той же 
центральной зоне равно 595 - - ліс, гдес 
-- плотность снега; влияние уровня грунто- 
вьшх вод на вариации До учитьваєт формула 
59 є дп Їсй ,гдей - - амплитуда колебания уров- 
ня вод; и, наконец,техногенное изменение рель- 
ефа -- соотношение 50 - у ЛП / | ,где АЛ -- 
превьшение между пунктами наблюдений за 
1 год, І. - - расстояние между пунктами. 

Необоснованное упрощение аналитических 
моделей средвбі с целью снижения неоднознач- 
ности истолкования во многих случаях явля- 
ется причиной неверньгх результатов вьгчис - 
лений геометрии источников, вертикального и 
латерального распределения плотностнь(х не- 
однородностей, особенно в случаях, когда сре- 
давокруг аномального источника далека от пред- 
положений об однородности. Надежное коли- 
чественное истолкование динамики масс, на- 
пример в мониторинге месторождений угле- 
водородов, получают при известной геомет- 
рий тяготеющих тел (по данньїм сейсморазвед- 
ки) и комплексной интерпретации поля силь 
тяжести и деформаций рельефа. 

Требование максимального учета имеющей- 
ся априорной информации о среде предлага- 
ем осуществлять двумя путями: 

1) путем построения надлежащих модель- 
ньх представлений (звездньге области извест- 
ной плотности внутри компактньх множеств 
в банаховом пространстве данньх); 

2) добавлением в функциональ типа невяз- 
ки в регуляризирующих алгоритмах некоторьх 
стабилизаторов дифференциального вида, соб- 
ственнье функции которьх совпадают о соб- 
ственньми функциями исходньх операторов. 

Обоснование п. 1 осуществлено в работе 
|Дубовенко, 2002), а п.2 -- в работе |Черная, 
1999). Решение конкретньгх обратньтх задач ме- 
тодом регуляризации осуществляєтся по алго- 
ритмам, изложенньх в работе |Тихонов и др., 
1983| и подобньх им -- на основе разработок 
|Дубовенко, 2002; Старостенко и др., 2004). 

О применениий мониторинга. Повторя- 
емье измерения гравитационного поля при- 
меняются в геодезиий, сейсмологий и вулка- 
нологии (комплексное определение предвест- 
ников землетрясений и извержений), для мо- 
ниторинга флюидной динамики водньх бас- 
сейнов (рис. 4), месторождений и подземнь(х 
хранилищ углеводородов в процессе их зкс- 
плуатации. При коррекции методики резуль- 
татьт мониторинга можно распространить и 
на сферу применения КГІВ, указанную в ра- 
боте |Собакарь, 19721. 
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При сравнениий данньх гравимониторин - 
га месторождения и его аналитической моде- 
ли обнаруживаєемое отклонение наблюденно- 
го поля от модели добьічи является сигналом 
к детальному изучению месторождения мето- 
дами ГИС и сейсмометрии. Из решения об- 
ратной задачи гравиметрии по данньгм мони- 
торинга можно восстановить общую картину 
вариаций плотности (снижение ДЛ? обознача- 
ет уменьшениєе обьема углеводородов вследст- 
вие их добьти и опускание газонефтяного кон- 
такта, увеличение Де -- поднятие уровня пла- 
стовьгх вод), но не абсолютнь6е значения плот- 
ностей. 

Гравитационньй мониторинг миграции СО» 
в месторождениях углеводородов основан на 
снижений их обьемной плотности при нагне- 
таний СО» в месторождениє. Так, согласно |Саз- 
регікоуа, Ноует5іеп, 2008|, ожнокилометровьій 
фронт флюмда СО, внутри 20-метрового пла- 
ста соли на глубине 1,9 км в условиях 30 9; 
насьщения СО» и 70 9; насьщения раствором 
соли генерируєт на поверхности гравианома - 
лию 10 мГал. Решение прямой пространствен- 
ной задачи с помощью |Старостенко и др. 2004| 
позволяєт по данньм микрогравиметрического 
мониторинга в процессе нагнетания жидкос- 
тей в подземнье хранилища вьіявить их рас- 
пределение и характер миграции по распре- 
делению плотностей. 





Глубина 
бассейна, м 











л9, мГал 





лу, мГал 








Рис. 4. Разрез бассейна с речньшми уровнями водьт 
и соответствующие вариации силь тяжести. 
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Большие перспективьт имеет решение со- 
пряженной обратной задачи -- истолкованиє 
временньх вариаций гравианомалий, возни- 
кающих вследствие изменения положения во- 
донефтяного контакта либо уровня пластовьіх 
вод в недрах либо скважинах. Его можно ис- 
пользовать как недорогой способ гравитаци- 
онного мониторинга подземнь(іх зкосистем ме- 
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